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, This paper discusses the elastic behavior of polyelectrolyte networks. The deformation behavior 

^\ ' of single polyelectrolyte chains is discussed. It is shown that a strong coupling between interactions 

0^ I and chain elasticity exists. The theory of the complete crosslinked networks shows that the Flory - 

Rehner - Hypothesis (FRH) does not hold. The modulus contains contributions from the classical 
^ ' rubber elasticity and from the electrostatic interactions. The equilibrium degree of swelling is 

0^ , estimated by the assumption of a c*-network. 
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^ ■ I. INTRODUCTION 



' The elastic behavior of polyelectrolyte networks is not understood. The main reason is the lack of knowledge in 
the interplay between the elastic part i^^i of the free energy and the part form the interactions i^int- The simplest 
approximation to overcome this problem is the application of the Flory - Rehner - Hypothesis (FRH) which suggests 
(— I , that both parts of the free energy can be added . Therefore the thermodynamic behavior of the network is then 
Q ' determined by the minimum of the total free energy F = F^x + Fi„t. This procedure must be wrong in general as 
O " pointed out by Ball and Edwards E|. The first reason is purely energetic. If the partition function is formulated by 
using the Edwards technique (see \M for a review), it must be expected that a cross term between both contributions 
\ exists. On the other hand it was shown recently that for realistic networks with fractal structural units, "left overs" 
^ . from the network formation processes yield immediately a violation of the FR - hypothesis Q. These results have 
' been shown by analytical considerations and by numerical simulations. 

Clearly the violation of the FRH was here shown to have purely structural reasons and it can be argued that 
for perfect networks the FRH should hold, as it was believed to be confirmed by many experimental work ||^,^ on 
^ I neutral networks and conventional rubbers upon welling experiments. Indeed it is difficult to say anything on the 

' validity of the FRH in swelling of conventional rubbers. Real rubber samples cannot be characterized as they should 

^ be and only indirect measurements such as neutron scattering [Q can give some hints, but the conclusions in favor 
"Y"; ■ or against the FRH seem to be vague. For a clear statement on the validity of the FRH, classical rubbers suffer 
d ' from the fact that the interactions in these systems are weak. The interactions are "only" excluded volume forces 
^ I and the corresponding potential is very short ranged. The short ranged excluded volume potential means that Flory 
I Huggins type, or Random Phase Approximations (see for example [^,|ll|) work very well. However, the excluded 
' ^ , volume interactions are weak and thus their effect on the elastic properties is indeed difficult to study with swelling 
experiments only. 

Polyelectrolyte networks do not suffer from the lack of weak interactions. Moreover the range of the interaction 
can be changed during experiments. Very often the interaction between chain segments are described by a Debye - 
Hiickel potential of the form V{v) cx l/rexp(— Kr), where r is the spatial distance between two segments and n the 
Debye screening parameter, which depends mainly on the total ionization / [p|-|Tot, i.e., oc /. The latter means 
$^ that the total ionization can be changed easily by adding and removing salt to / from the solution. The range and 
5^ , the effect of the interactions can be changed easily. Large ionization / corresponds a high concentration of salt and 
thus to large values of k. In this case the interaction potential can be safely modeled by a short ranged form, in some 
sense with an excluded volume type potential with a certain value determined by the Bjerrum length and k p^ . The 
more interesting case is therefore the regime of small values of the Debye screening parameter corresponding to long 
range potentials. This is the regime discussed below in more detail. 
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We will show below systematically that the FRH is not applicable to the case of polyelectrolyte networks. To do so 
we first derivate the Flory-Rehner Hypothesis. Then to make the paper self contained we start with the computation 
of the deformation of a sin gle polyelectrolyte chain, though some progress has already been achieved by the 
application of blob-pictures . Indeed these ideas have been also applied to networks . The problem is here is 
that the such obtained results are based on single chain approaches. Therefore these results are restricted to cases 
where either the meshes are largely stretched, or the electrostatic interactions are screened. These two cases are 
reflected on the ratio of the screening length and the chain size, i.e., 1/k much larger than the mean mesh size radius 
the network or vice versa. In the latter the network behaves almost as a classical swollen network, whereas the first 
case corresponds to the strong polyelectrolyte regime jl^ . In the present paper we will also present a theory beyond 
these two extreme regimes. Moreover we will leave in the following the single mesh approximation and study a network 
composed by many meshes. Then we cannot rely on the extreme cases k large or small.. 

Our results will show that there is a strong interplay between the elastic force and the interactions. The force 
extension relationship is in contrast to neutral Gaussian chains not only dependent on the chain length and temperature 
I psf but also a strong effect on the single chain modulus from the electrostatic interaction appears . This method 
will be generalized to polyelectrolyte networks jl^]. We show below that the elastic modulus is a combination of 
both contributions, i.e., the classical elastic part, determined by crosslinks only | [T^ and a new part coming from 
the interactions. Thus the modulus is determined strongly by the strength and the range of the interaction. Note 
that this is not the case for neutral networks. All computations for excluded volume networks so far show, that the 
modulus is determined by the number of crosslinks only [p^ , p^ . In the third main part of the paper we use the single 
chain results to study the swelling behavior of polyelectrolyte networks. Before we proceed we present a mathematical 
formulation of the origin of the FRH to show what has to be done to go beyond the classical theories. 



II. DERIVATION AND LIMITS OF THE FLORY-REHNER-HYPOTHESIS 



We consider a polyelectrolyte network in solvent containing A^i indistinguishable counter-ions. The partition function 
for such a network consisting of one very long single chain with a fixed crosslink-configuration S in terms of a path 
integral reads: 



Z(S) = 1 Vv{s) I^...J^ eM-m n S[r{s,) - r(s,)]. 



(1) 



here denotes that the i-th and the j-th monomer form a crosslink, r(s) represents the chain conformation in 

three dimensions as a function of the contour variable s and R; is the position of the l-th counterion. Vq denotes the 
integration-volume. As we assume that the network is weakly charged the underlying Hamiltonian H is given by an 
Edwards Hamiltonian for flexible chains. Considering the counterions explicitly the Hamiltonian H is: 
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where / is the Kuhn length, A'o the total number of monomers of the crosslinked chain, /3 is {kBT)~^, where fee is 
the Boltzmann constant and T denotes the absolute temperature, eo is the dielectric constant and er the relative 
dielectric constant, z is the monomer charge in units of e and is charge of the i-th counterion, which is assumed to 
be independent of i. Ri denotes the position of the i-th counterion. 

In order to use collective coordinates, we introduce the densities pp of the polymers and pi for the counterions. The 
constrained partition function for fixed densities pp and pi reads (see Appendix A): 



ZiPp,PuS) = exp 



X exp 
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Vpp{k)j pp{k)pp{-k) 
VCo(k) ' ^i(k))pi(k)Pi(-k)-5]ypi(k)pp(k)pi(-k) 



(3) 



with Vpp(k) = 5z^/k^. Vpi and Va are defined analogously. S'o(k) and Co(k) are defined as: 
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5o(k)=(/ ds ^g' ^Mr(s)-r{s')) \ 



and 



where < . . . >o means the expectation value with respect to Hq 



Co(k) = ( ^e^'^^^-^'A (5) 

' 



2P Jo \ds 
Eq. (P) defines an effective Hamiltonian Hcg: 



/5^o = 4 r'ds f In \ II S[r{s.)-r{s,)] (6) 



/3i?cff(Pp, P^^^IT. + ^PP^'^)) ^p(k)Pp(-k) + i + Mi(k)) Pi(k)pi(-k) 

+ ^Fpi(k)pp(k)pi(-k) (7) 

k 

From the effective Hamiltonian we separate the k = term, i.e.: 



/377eff = /3i/°ff(Pp(0),P^(0)) +^/3i/^ff(pp(k),p,(k)) (8) 

A:>0 

with 

/3i7°ff(pp(0),p,(0)) - I + V(0)) Pp(0)Pp(0) + i + Mi(0)) a(0)a(0) 

+ Fpi(0)pp(0)A(0) (9) 

For pure Coulomb potentials in infinite large systems the values for Vpp(O), T/pi(0) and Vn{0) appear to be divergent. 
In realistic systems the chains (and meshsizes) are finite. The lowest value for the wave vector is then of the order of 
|k-^| ^ \/{1'^Nq) which acts as a cut off for small wave vectors. Thus divergences are prevented. 

If fluctuations are neglected the only contribution to the partition function stems from . In this case the Free 
Energy Fq reads: 

Fo = -ln(Zo(pp(0),pi(0))) (10) 
= \ + + Pm ln(p,(0)) + il^pp(0)p^(0) + Iv^mPiiO) + l^pi(0)pp(0)pi(0) 



In deformed systems the first term of the right-hand-side of Eq. (|1C|), which is written in brackets, contains the elastic 
Free Energy. The second term represents the translational entropy of the counterions and the last three terms are 
interaction parts. This equation indicates, that neglecting fluctuations the Free Energy is just a sum of an elastic, 
an entropic and an interaction part F = F^i — TS'countor + -Pint- This reflects the Flory-Rehner Hypothesis. It shows 
up only on the limit k = 0. Below we will show that this is not all all sufficient. In general fluctuations cannot be 
neglected. In the present network case fluctuations have a strong effect: they violate the Flory-Rehner Hypothesis 
and renormalize the Gaussian elasticity. If the Free Energy is calculated without neglecting the fluctuations one 
has to integrate over all jOp(k) and pi(k) by performing the path integrals. Eq. (|^) shows that in this case there 
will be a strong coupling between the elastic and the interaction part of the Free Energy. Even in a Debye-Hiickel 
approximation, i.e. pp(/c > 0) = this coupling cannot be avoided. In the following sections we show, that in 
polyelectrolyte networks there is a coupling between the elastic and the interaction part and thus the Flory-Rehner 
Hypothesis is no longer applicable. 



3 



III. SINGLE POLYELECTROLYTE CHAIN 



To approach the problem of elasticity in polyelectrolyte networks on a systematic basis we recall the results of the 
single chain deformation. This section will introduce the mathematical method and - more important - it will be 
shown that the force - extension relationship is fundamentally changed from what is known in Gaussian and excluded 
volume chains. The deformation of the end - to end - distance of a Gaussian chain by an external force / is simply 
described by 

/ = fcBT^ (11) 

where T is the temperature, I the Kuhn length, N the chain length and r the distance of the chain ends. This result 
can be easily derived from the Gaussian statistics of the ideal chain Already when excluded volume forces are 
taken into account eq.(^l]) is changed in a fundamental way. Excluded volume force change the statistical behavior 
of the chain, i.e., the size if the chain is given hy R~IN^/^. Thus the single chain is swollen, and another stretching 
law is expected. This problem was solved by Pincus The result is briefly quoted. The problem is no longer 

isotropic and the force has two components, i.e., parallel to the stretching direction and perpendicular to the stretching 
direction. At low deformation a Hookian deformation behavior is expected, whereas at larger deformations the self 
avoiding statistics comes into play. In the low deformation regime the result is 

l(r||)|~(R')/ (12) 

and in the limit / ^ oo 

|(r||)|^7V/i/^-i (13) 
for the elongation in parallel direction to the force f . For the perpendicular direction it is found (/ S> 1) 

(ri) ^ Ar/i/-2 (14) 

A similar situation will occur in the case of polyelectrolyte chains. The coupling between elasticity and interactions 
will become much stronger. The interactions are very strong and pre-stretch the chain to a large extend. Thus it is 
instructive to study the elastic behavior of single polyelectrolytes. 



A. Model 



Let us first introduce the standard model which is employed here. Since we restrict ourselves to flexible chains 
which are weakly charged the Edwards model is the appropriate tool, like already given in section II. Here we 
consider a Debye-Hiickel chain, thus neglecting fluctuations of the counterions: 

3 /■^'' . /dr\' . /-^^ . „dr 



Wr;f]^^y^ d,s ^^j '^'d. 
hz^ /■^« /-^"^^ , exp{-^|r(s)-r(,sO|} 

where b = e^/47reoerfcBT' is the Bjerrum-length, e is the electrical charge of an electron, eo is the dielectric constant 
and er the relative dielectric constant, f is the external force and k^^ denotes the Debye-Hiickel screening length. 
The correlation function can be calculated in terms of a path integral [|l]-|23j as follows: 

<.r(Afo)=r 

G(r,iVo;f)= / I?r(s) exp {-/3i/E[r; f]} . (16) 

Jr(0)=O 

Its Fourier transform is defined by 

G(k, No; i) = J d^r exp{-ikr}G (r. No; f ) . (17) 

The averages of the force-size relationship (R-^) and (R'^), where R|| denotes the parallel component with respect to 
f and R_L is the corresponding perpendicular part, are then readily calculated by the general formulae 
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and 
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(Ri)[f] = - 



d^dki G(k,7Vo;f) 



G(k,7Vo;f) 



k=0 



Et=idVdkf G(k,7Vo;f) 



G(k,7Vo;f) 



(18) 



(19) 



k=0 



By analytic continuation of the Fourier space to the complex plane, the correlation function G(k, iVp; f) can also be 
written as the zero- force correlation function G(k — i^f, -/Vq; f — 0). Substitution of Eq. |l^ into Eq. |l^ yields: 



G(k,7Vo;f)= J d^r exp{-ikr} 
For constant force f Eq. (EH) can be rewritten as: 



r{No)=r 



r(0)=0 



Vv{s) exp{-/3i7E[r;f]} 



(20) 



G(k,iVo;f)= J d^r exp{-i(k - i/3f )r} 
= G(k~ip{,No;{ ^ 0) 



r(Aro)=r 



r(0)=0 



Vr{s) exp{-/3ifE[r;f = 0]} 



(21) 



Consequently, to get results for (R.|) and (R'^), we only have to calculate G(k, A'o) and continue the first argument 
of G to the complex plane. 



B. Calculation of the Correlation Function 



For technical reasons it is convenient to introduce a field theory, which enables to carry out the variational technique. 
Thus the field theoretical Hamiltonian reads |16|: 
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V'(ki)V'(k2)(7(ki + k2)(5(ki + k2 + ka + U)i'{^z)^{\^4.) 
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where J^. is an abbreviation for / (P'k/ (27r)^ and U denotes the Debye-Hiickel potential in units oi (3 ^. In the Fourier 
space G(k, /io), which is the Laplace transform of G(k, A^o) with respect to A^o, can be written exactly as ||2^ : 



G(k,/xo)= (^Mo + ^fc' + S(k)^ 
where E(k) denotes the proper self energy. We now define a trial-Hamiltonian Ti: 

where ^(k, /io) is an approximate correlation function with an approximate proper self-energy M(k): 

g{k,fio)= (^^io+^^k^ + Mik] 

In this notation the well-known Feynman inequality is given by: 

F<T+{H-H)n 

where 

/P^...exp{~/37^} 
(. . .)-H = hm 

n^o J x)^ exp{-/3H} 



(23) 



(24) 
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is the mean-value and T the free energy with respect to W. The right hand side of the inequahty (|2^) has to be 
minimized with respect to M.^ which is our variational function. The general minimization condition reads 

\T +{H-n)n) = ^ (28) 



(5A^(q) 

where (5/5A4(q) denotes the functional derivative with respect to A^(q). This condition ( [2^ ) is equivalent to ||l6| 

-^(q) = T 



k MO + f fc2 + A^(k) 



6 ' 

4^^^^' / . , V , (29) 

"6' 



k {k^ + (q + k)2) (;xo + f fc2 + Al(k)) 



This is a non- linear integral equation for Al(q), which in the following has to be solved approximately, since the 
exact solution is unknown. At this point it should be stressed, that Eq. (|2^) represents the well known Hartree 
approximation. As our aim is to get a force-size relationship according to Eqs. (|l^) and (|l9|) the ansatz for the 
approximate proper self energy 

M{ci)^M (0) + ak2 -f ©(k'*) (30) 
is valid for small external forces /. The constant A^(0) is absorbed by the shifted chemical potential ^ — hq + A^(0). 



Consequently the only task is to calculate the coefhcient a self-consistently from Eq. (29). This is possible for small 



inverse screening length k, which means in the long ranged limit of the Debye-Hiickel potential. The result is: 



2bz^ r 

a 



For details see p^ ]. 



l + O 



kI 



(31) 



C. Results 

Let us describe the consequences briefly. Inserting the approximate result for the proper self-energy S(k) into 
Eq. ( ^ ) yields an explicit expression for the correlation function G(k, fi) with the shifted chemical potential fi as 
mentioned above 

G(k,M)=fM+?fc^ + |^fc^) ' (32) 

Now the conformational free energy of the chain under the influence of a force f can be calculated very easily by 
0(/i, f) — — lnG(0,M;f) = — lnG(— i/3f, /i). Using the well-known thermodynamic relationship 

we express fj, depending on its conjugate variable iV and the force f . Note that the variable N is not the bare number 
of monomers Nq since ^ is a. shifted chemical potential, but N is proportional to A'o, indeed it is easily seen that 
N < Nq. Considering only singular terms in k and neglecting terms of order this calculation yields: 

M=^+ ' +0{f) (34) 

Substituting Eq. (||) into Eq. (H) we get G(k, A^). According to Eqs. (|l|) and (|l9|) (Rjj) and (R2_) can be calculated 
from Eq. (|3^) . Expanding in a power series for small forces to second order and again considering only most singular 
terms for small k, (R|) becomes for /3//k < 1 
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The square root of the mean square elongation can be written according to Eq. (Bq) as 



;Rp[/]-(Rf)[0]~/ (36) 
which is a Hook-Hke law. 

Using the same approximation as mentioned above the root mean square end-to-end distance perpendicular to the 
force /, ^ (R-j_), decreases with / for /3//k < 1, which is contrary to a Gaussian chain In particular: 



Ri>[/] = ^^3- - /^^/^^J-J-^ + (37) 



For f — the perpendicular end-to-end distance squared becomes doubles exactly (R|)[0], since we have assumed 
rotational isotropy. 

The latter relation ship for zero forces does not agree with one of the two scaling limits, i.e., strong polyelectrolyte 
regime or excluded volume regime (large k). In the strong polyelectrolyte regime, very small k the chain size does 
not depend on the screening parameter p^ , whereas here it does. The reason is that we are here in the intermediate 
regime where 1/k is of the same order as the chain size. This is visible by a detailed study of the renormalization in 
the chain length from iVo to N. It is exactly this regime which we need to study the statistical mechanics of networks. 

IV. POLYELECTROLYTE NETWORKS 

In this section we calculate the modulus of a polyelectrolyte network with a random crosslink configuration. 

A. Model 

A network is formed out of a very long chain by the instantaneous introduction of a sufficiently large number of 
crosslinks in the liquid phase (see for example We restrict ourselves to a network, which consists of flexible, 

weakly charged strands. Consequently the Edwards model is an appropriate tool to describe the network 
Therefore the underlying Hamiltonian is the same as for the single chain Eq. (^5|). Again we assume Debye-Hiickel 
interaction between the monomers. For the introduction of crosslinks we choose the standard way suggested by 



Deam and Edwards |18|. We assume for mathematical convenience four functional crosslinks which join two arbitrary 
segments r{si) and r(sj) along the chain. Of course, the value for the free energy then depends on the specific choice 
of the pairs of monomers, but on macroscopic scale only the statistical average on any crosslink configuration is of 
importance. Nevertheless this requires non - Gibbsian statistical mechanics in the sense that the crosslink positions 
represent quenched degrees of freedom. 

To calculate therefore the free energy F of the network, we have to take the statistical average over all crosslink 
configurations S. This represents the fact, that -F is a self averaging quantity. 



FiNtot,N,) = -kBT J dSV{S)\iiZ{S). (38) 

Z{S) is the constrained partition function for a network with the crosslink configuration S (see section II), Nc is 
the number of crosslinks and 7^(S) is the crosslink distribution function. Since we assume that the crosslinks are 
instantaneously introduced in the liquid phase, V{S) is yielded by the constrained partition function of the liquid 
phase Z(°)(S) Eq. ([|). The crosslink distribution function 7^(8) is simply: 

In the following it appears to be reasonable to assume that the so chosen distribution function does not depend on 
the specific deformation of the network. Note that ^(S) differs generally from Z^°\S). Since we are interested in 
deformations of the network, Z(S) is the partition function of the deformed network. 

To calculate the free energy F (Eq. (^)) explicitly it is convenient to use of the so called replica trick [Q. Define 

FMouNc) = -A^bTIu J dSZ(")(S)Z'"(S), (40) 
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where m is the Rephca index. The free energy F, which is averaged over the disorder of the crosshnks, reads ||l 

dm 

As in the case of the single chain, the free energy F is calculated by making use of its relation to the corresponding 
distribution functions and Green functions of the corresponding propagator (see for the technical details). 



(41) 



B. Field-theoretical Formulation and Results 



As in the case of the single chain the external force acts on the ends of the chain and is reintroduced by the analytic 
continuation to the complex plane: 



G(k, ^0, 2c, f ) = G(k('') , k(i) - i(3f, . . . , k(") - if3{, ^, z^, f = 0) 



(42) 



Here G'(k, /io, Zc, f ) is the Greens-function, where k is the wave vector in the 3(m + l)-dimensional Fourier transformed 
replica space and k'*) is the three-dimensional wave vector of the i-th. replica, /io is the chemical potential of the 
monomers and Zc is the fugacity of the crosslinks. To calculate the Greens function G(k, /xq, Zc) we introduce a field 
theory. The field theoretical forceless Hamiltonian for the network in Fourier space reads Eol], |T6|: 



Him)] = 



|V(q)^(-q) + ^q''/'(q)'/'(-q) 



■0(qi)V'(q2)V'(q3)^(q4)'5(qi + q2 + qa + q4) 



qi.q2,q3.q4 



E 

fc=0 



qi,q2 



V^(qi)V^(q2)n'5(qf^+q^'^) 



q3,q4 



V^(q3)v^(q4)n'^(q3^ + q4^) 

l^k 



(k) 



(fc) 

qa 



qf) 



(43) 



where F^'"'^(g) is the Fourier transform of the Debye-Hiickel potential in the fc-th replica and an n-component vector 
field. The further procedure of the calculation is exactly analogous to the calculation concerning the single chain 
presented above. Again we introduce a Gaussian trial Hamiltonian with the correlation function: 



^(k, /io,Zc 



Mo + -fc^ + Al(k,zc) 





(44) 



Then we make use of Feynman's inequality for a variational principle. The variational parameter again is the approx- 
imate proper self energy Al(k, Zc). The only difference is, that now the proper self energy depends on a replicated 
Fourier vector k. For details see 17 1. Assuming a replica symmetry this procedure leads to the following self-consistent 
equation for M(k, Zc) 



M(k,Ze) = 



2Trbz^n{m + l)Vo" 



1 



Airbz'^r 



q M + |q^ 

1 



M(q,Ze) 



1 



q(o) (q(0) +k("))2 + k2 /i+ |(q(0),k(i-™))2 + M(q(0),k(i-™)) 
1 1 
(1) (q(i) +k(i))2-f-K2^+ |(k(0),q(i),k(2-™))2 -|-M(k("),q(i),k(2-™)) 
1 



(45) 



q M + 



M(q,2c) 



where k*^' '™-' = (k^*\ . . . , k*^™^). This equation has to be solved approximately. In analogy to the calculations on the 
single chain pa we restrict ourselves to small external forces applied on the ends of the crosslinked chain. Therefore 
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we make the same ansatz for the proper-self energy as in the case of the single chain. The only difference is that the 
proper self energy in this paper is a function depending on the replica-space wave vector q: 



Af(q) = ao + aiq2 + 0(q^) 



(46) 



In the limit of a weakly charged network and in the long ranged limit of the Debye-Hiickel potential the result of the 
self-consistency calculation according to Eq. (Ea) for gq and ai is (for details see [O]): 



ai 



2bz^ 



l + O 



kI 



and 



ao 



/i \ TT 2 9/i7r 

K f2K^Zcm K'^Zcjm k^Zc ln(fi;//y^)rn k^Z; 



(47) 



(48) 



37r2 



47^2 



27r2 



67r2 



Therefore the approximate Greens function reads: 

^(k,/i,Zc) = 



/i + ao + f-k2 + aik2 



(49) 



From Eq. (^) the grand canonical partition function in replica space under the influence of an external constant 
force / acting on the ends of the chain can be calculated: 



Erni^i, Ze, f ) = ^(k(°) , k^^) - Z/3f , . . . , k^") - l/3f , fi, Ze) Ifc^Q 

1 



fj, + ao - ^m/32f2 — aim/32f2 



(50) 



Here we reintroduced the force / according to Eq. (|4^). The monomer chemical potential and the fugacity of cross- 
links which parameterize the grand canonical partition function in Eq. (^), should be expressed in terms of the 
parameters iVtot and Nc. According to Panyukov and Rabin the expression Fm{Ntot, N^) can be calculated by 
the method of steepest descent in the thermodynamic limit iVtot, — oo: 



FmiNtot,N^)/kBT = - lnS(/i, z^) - iVtot/i + A^c InZe 



(51) 



Consequently the fugacity z^ of cross-links and the chemical potential fj, of monomers can be obtained by minimizing 
the right-hand side of Eq. (pl|): 



and 



Nr. = 



din ^,n^o(/i,Zc) 

d\n S„^o(/i,Zc) 
91nZc 



(52) 



(53) 



in the limit of a vanishing replica index m. In the weakly charged limit of the network and in the long-ranged limit 
of the Debye-Hiickel potential the result for yu and Zc is: 



and 



fi = N 



67r2iV2 37r(4 - 7r)&z2 



Moreover it can be shown [n9|, that the conformational free energy of the network is given by: 



(54) 



(55) 
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dm 



(56) 

m— 



If the free energy F is known the force-size relationship is the simply calculated by the derivative of F with respect 
to the external force /: 

(i?) = ^ "'■^^ (57) 

After having inserted the results for Zc and /i in the force size relationship, expanded for small charges z neglecting 
terms of order and higher, since only weakly charged networks are stable, and considered only terms of leading 
order with respect to small k the force-size relationship reads: 

This is again the result for the small deformation regime. The result describes a Hookian law for the force extension 
relation and defines the elastic modulus of the polyelectrolyte network. Note further that this result is valid for small 
forces, i.e. /3//k < 1. Therefore the modulus for the small screening and the low deformation regime of the network 
reads: 

iWh^ ^ 6iV 



(59) 



The modulus depends on the density of the crosslink and on the Debye screening parameter. Thus both contributions 
enter in a significant way. Most important is that part of the modulus stemming from the interactions depends on 
the crosslink density iV squared. This now indicates that in the case of strongly interacting soft materials the FRH 
is no longer applicable. The modulus computed here is a strong combination of both contributions. Moreover we 
stress that the result is beyond perturbation theory and RPA - type approximations in contrast to earlier suggestions 
published previously 123]. The author in this paper suggests a modulus which is not changed by the interactions. 



V. SWELLING BEHAVIOR 



Let us finally discuss the consequences of the above results to the swelling behavior of a simple polyelectrolyte 
network. The general theory of swelling for these systems is not simple and will be studied in a subsequent publication. 
Nevertheless we have a good basis to study the general swelling behavior of polyelectrolyte networks. To do so we 
assume that the equilibrium swelling is given by the c* - network jsj . The consequence of this assumption is that the 
single chain behavior can be used for the determination of the degree of swelling. Thus we may restrict ourselves on 
one mesh of the network, and we can use the single chain behavior studied in the first section of the paper. Of course, 
we have to neglect complications by the presence of entanglements and topological restrictions. Therefore we can 
apply the following results only for mesh sizes below typical entanglement lengths. The strong electrostatic repulsion, 
however, does not allow for a large degree of topological restrictions, which makes the assumption of no entanglements 
reasonable. To estimate the equilibrium degree of swelling, we use an ideal c* - network whose monodisperse meshes 
consist of chains with N monomers. These chains are connected by the process of endlinking, where the functionality 
of an end of the chain should be six in three dimensions. Thus in the ideal totally stretched state of the network the 
end- links are on the lattice points of a simple cubic lattice with the lattice constant Nl. Figure 1 shows such a swollen 
network in two dimensions. 
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FIG. 1. Swollen network in two dimensions 



For long chains this ideal network is of course not very realistic. On the one hand entanglements occur, which 
disturb the ideal structure drastically. On the other hand two points of the network could be connected, which would 
be in the ideal, swollen state far away from each other. If the network, however, is synthesized out of short chains, 
that means N is small, then these defects should occur with a much smaller probability. Consequently in this section 
we have to restrict ourselves to networks, which are made of very short chains. 

In the dry state without having contact to a solvent the considered network is a dense system of neutral polymers, 
because counterions are bound to the charges on the polymers. In solvent some counterions dissociate and as a result 
we get a screened Coulomb interaction between the charges, which swells the network. To describe the swelling of a 
network correctly, one has to know the complete free energy of all the states between unswollen and swollen. This 
means that we have to have a theory that is valid for the network from the unswollen to the swollen state. As far 
as we know such a description does not exist by now, because theoretical considerations are only made in special 
physical limits. Therefore in this section we use the idea of the c* - network: we compare the volume of the dry and 
the swollen state of the above mentioned model-network. 

The dry state is a dense system of neutral polymers, which means that the remaining excluded volume interaction 
is screened [0. Consequently the average end-to-end distance of such a chain is: 

(Rl) cx PN (60) 

We assume that the considered model network is formed in an equilibrium state, where the mean end-to-end distances 
of the chains obey Eq. (pQ). In the swollen state the averaged end-to-end distances of the chains connecting two 



neighboring endlink-points in the network can be calculated from Eqs. ( pSD and (37). For vanishing external forces /: 

(X (61) 

Thus for our model network we can introduce a parameter of swelling Q defined os Q = Y jV^^ where Y is the volume 
of the swollen network and Vq is the volume of the unswollen network: 

Here N — Nc/Ntot is the number of endlinks Nc divided by the total number of monomers iVtot, i. e. the endlink- 
density. In our case N — 1/N is the averaged number of monomers per strand. 

Here we see clearly that in contrast to neutral networks the degree of swelling Q does not only depend on the 
crosslink density, but also on the range of the interactions via k and therefore the ionization. A typical dependence 
on the different quantities is shown in figure 2. 
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FIG. 2. The degree of swelling as a function of crosslink density (dashed line) and Debye screening parameter (bold-faced 
line) . 

VI. DISCUSSION 

In this paper we developed different stages of the problem of the elastic behavior of polyelectrolyte networks. This 
problem required a new approach and the classical theories for neutral networks cannot be used, so far, except in the 
limit of large salt concentration, i.e., k oo. In this regime the interaction potential between the segments is very 
short ranged and the network behaves roughly as an " effective excluded volume network" . Naturally this regime is 
not very interesting, since no real new insight can be expected psf . On the other hand, the small k regime is very 
important and novel theories can be formulated In this regime wc found the strong invalidity of the Flory - 

Rehner - hypothesis, which states that different parts of the free energy can be added. The elastic modulus of the 
network contains both: the entropic part and the interactions. Both seem to show up as two springs in a series. The 
equilibrium swelling is determined by interactions and crosslink density. 
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APPENDIX A: 

Starting from Eq. (|^) we introduce the densities Pp(k) of the polymers and pi(k) of the counterions. The result is: 



Z(pp,pi,S) 




with b = / 4:7760 f-rkBT the Bjerrum-length. According to the (5-distributions in Eq. ( [Al| ) the partition function can 
be rewritten as: 
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X exp 



-J E ^pp(k)Pp(k)pp(-k) - i ^ l^i(k)pi(k)pi(-k) - J2 l^pi(k)pp(k)pi(-k) 



where Vpp(k) = 6z^/k^. Vpi and Vii are defined analogously. The (5- functions in Eq. (A2) can be rewritten by 
a (h- and (A'-field: 



dR 



l?0(k) exp 



1 



*0(k) Up(k) - dse*--(«) 



X / V<j)'(k)exp 

3 /-A^o 



X exp 



X exp 



2P 



ds 



^) -ln|n^[r(^.)-r( 



i ^ ypp(k)pp(k)pp(-k) - i ^ ^^i(k)pi(k)pi(-k) - ^ ypi(k)pp(k)pi(-k) 



k k 

Within a Gaussian approximation the partition function reads: 



,»0(k)pp(k) 



2 JO 

Ni 



No nNo ^ 

ds / ds'e'''('-(^)-'-(^'» 



y P,?i'(k)e*'^''''''('^) ( l-i(/)'(k)^e 



X exp 



X exp 



' ( ^ - In [] 5[r(s,) - r(s 



1 
2 



2/2 7p " \dsj 

-\ ^ ypp(k)pp(k)pp(-k) - 1 ^ l^i(k)Pi(k)pi(-k) - ^ l/pi(k)pp(k)pi(-k) 



Denote by S'o(k) and Co(k): 



and 



Soik) = ( / ds ds'e*('-(^)''-(^» 
\Jo Jo , 



Co(k) = ( e''^^^^ -^') 



where < . . . >o means the expectation value with respect to iJo with: 
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Line two and three of Eq. (A4) are now replaced by the mean values of the terms concerning _ffo- The result is 
considered as the first terms of a series expansion of an exponential function. Consequently the exponential function 
is reintroduced. After having integrated over and (j)' the partition function reads : 



Z(pp,pi,S) = exp 



X exp 




Co(k) 



5o(k) 



1 



1 



+ V^pp(k) pp(k)pp(-k) 



+ Mi(k) Pi(k)A(-k) - ^ ypi(k)pp(k)pi(-k) 



k 



(A8) 
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